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"^vq Abstract 

We assess the energy limit up to which Heavy Baryon Chiral Perturbation Theory can be accurately applied to the process of neutral 
^vq pion photoproduction from the proton by analyzing the latest data from the A2 and CB-TAPS collaborations at Mainz. We find that, 
within the current experimental status, the theory works up to ^170 MeV. Above this energy the data call for further improvement 
in the theory such as the explicit inclusion of the A(1232). We also find that data and multipoles can be well described up to ^185 
|^ MeV with Taylor expansions in the partial waves up to first order in pion energy. 

CO Keywords: Chiral perturbation theory, effective field theory, pion photoproduction, heavy baryon 
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1. Introduction 

Chiral Perturbation Theory (ChPT) is an effective field the- 
ory (EFT) of Quantum Chromodynamics (QCD) in the low- 
energy domain where quarks and gluons are confined into hadrons 
and conventional perturbation theory cannot be directly applied. 
Due to the spontaneous breaking of chiral symmetry in QCD 
the ir meson appears as a pseudoscalar Nambu-Goldstone bo- 
son {If] becoming the carrier of the nucleon-nucleon interac- 
tion. However, when fully relativistic spin- 1/2 matter fields 
(i.e. nucleon) are introduced in the theory the exact one-to- 
one correspondence between the loop expansion and the ex- 
pansion in small momenta and quark masses is spoiled |01. 
This is due to the fact that the nucleon mass M does not van- 
ish in the chiral limit. A consistent power counting scheme 
known as Heavy Baryon Chiral Perturbation Theory (HBChPT) 
3tl overcomes this difficulty considering the baryons as heavy 
(static) sources. For irN scattering and photo-pion production 
HBChPT has been very successful at describing experimental 
data in the near threshold region Q3|, |5|] . In this Letter we address 
the question of how well it works for the latest and most accu- 
rate 7p — > tt°p data to date [4] and to provide the first upper 
limit of accuracy for HBChPT in describing this reaction. The 
recently completed Mainz data is for the differential cross sec- 
tions da I dfl and linear polarized photon asymmetries £ for the 
7P — > ir°p reaction taken from threshold through the A(1232) 
region. This was performed with a tagged photon beam with 
energy bins of 2.4 MeV. We also determined the low-energy 
constants (LECs) to see if they are actually constant as the pho- 
ton energy is increased. The quality of the HBChPT fits -% 2 
per degree of freedom (x 2 /dof)- are also compared to a sim- 
ple empirical benchmark fit, a Taylor expansion of the partial 



waves. The data in J3l are at present the highest quality data 
available on neutral pion photoproduction from the proton and 
allow to perform analyses not possible before on the limits of 
chiral symmetry like the one presented in this Letter. 

2. Theoretical Framework 

Due to the symmetry breaking, the S-wave amplitude for 
the jp —> ir°p reaction is small in the threshold region, — van- 
ishing in the chiral limit fl. Additionally, the P-wave ampli- 
tude is large and leads to the A(1232) resonance at intermediate 
energies J6[]. Hence, for the -/p —> ir°p reaction the S- and P- 
wave contributions are comparable even very close to threshold 
J3l and even D waves have an important early contribution due 
to the weakness of the S wave [8]. The differential cross section 
and photon asymmetry can be written in terms of electromag- 
netic responses 
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where Wt and Ws are the electromagnetic responses, 9 is the 
center of mass scattering angle, fc 7 the center of mass photon 
energy, q the pion momentum in the center of mass, and s the 
squared invariant mass. The responses Wt and Ws are defined 
in term of the electromagnetic multipoles: 

W T = T {s) + T x (a) Pi (6) + T 2 (s) V 2 {0) + ... (3) 

W s = So (s) + Si (s) V 1 (9) + ... (4) 
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where Pj (9) are the Legendre polynomials in terms of cos 0, 
the dots stand for negligible corrections, and 

T n (s) = ^2Re{M*(s)T^M j (s)} (5) 

ij 

S n (s) = Yl Re i M * M S n M J ( s ) } ( 6 ) 

ij 

where M 3 (s) = E +, E 1+ , E 2+ , E 2 _, M 1+ , Mi-, M 2+ , 
M 2 _ . The coefficients and S% can be found in Appendix 
Ain@]. 

The partial waves (electromagnetic multipoles) are not ob- 
servables and have to be extracted from the experimental data 
within a theoretical framework. In this Letter we employ three 
approaches to describe S and P waves that we present in forth- 
coming paragraphs: Section|2j]HBChPT JnJGjl]; Sectionl2~2l 



Unitary HBChPT (U-HBChPT); and Section^ Empirical. In 
all cases D waves are incorporated using the customary Born 
terms. Higher partial waves can be safely dismissed in this en- 
ergy region H. The conventions employed in this Letter and 
further information on the structure of the observables in terms 
of the electromagnetic multipoles can be found in |9J]. 

2.1. HBChPT 

The explicit formulae for the S and P multipoles to one loop 
and up to 0(q i ) can be found in |T(1 U|. Due to the order- 
by-order renormalization process six LECs appear: a\ and a 2 
associated with the Eq + counter-term: 



Po+ = eaiwm^o + ea 2 U} 3 , 



(7) 



where oj is the pion energy in the center-of-mass; b p associ- 
ated with the P3 = 2Mi + + Mi_ multipole together with 
£1 and £2 associated with Pi = 3Pi+ + M\+ — Mi- and 
P 2 = 3£i+ — Mi + +Mi_, respectively. The C4 LEC associated 
with Pi, P 2 , and P3 has been taken from lfl2ll where it was de- 
termined from pion-nucleon scattering inside the Mandelstam 
triangle. Some other parameters appear in the calculation, but 
these are fixed. The full list is: the pion-nucleon coupling con- 
stant g^N = 13.1; the weak pion decay constant f n = 92.42 
MeV, together with the anomalous magnetic moments of the 
proton and neutron, the nucleon axial charge (which we fix 
using the Goldberger-Trieman relation cja = QttN fn /M); and 
the masses of the particles. The pair (ai, a 2 ) LECs are highly 
correlated, r(ai ,a 2 )= —0.99 ifsL 1 1 ill , and it is more convenient 
to use the pair of LECs (a + = a\ + a 2 , a_ = ai — a 2 ), where 
a + is the leading order for the counter-term close to threshold 
{uj ~ m^o) Jit]. Henceforth, five LECs are fitted to the data 
under this approach: a + , a_, £1, £2, and b p . 

2.2. U-HBChPT 

From general principles such as time reversal invariance and 
unitarity the S wave can be written as the combination of a 
smooth part and a cusp part [HBO 
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where 8q is the 7r°p phase shift (which is very small), yfs is the 

invariant mass, \J sj^™" 1 the invariant mass at the 7r + n thresh- 
old, q + is the ir + center-of-mass momentum, Aq is Eq + in 
the absence of the charge exchange re-scattering (smooth part), 
and j3 = ReEo+ (7P — > n + n) x a (7r + n — > ir°p) parameter- 
izes the magnitude of the unitary cusp and can be calculated 
1 13] on the basis of unitarity. Equation dHJ takes the static 
isospin breaking (mass differences) as well as irN scattering 
to all orders into account. In the electromagnetic sector it in- 
cludes up to first order in the fine structure constant a. The 
7r + center-of-mass momentum, q w +, is real above and imagi- 
nary below the tt + threshold; this is a unitary cusp whose mag- 
nitude is parametrized by j3 which can be calculated 111 311 on 
the basis of unitarity and taking into account isospin break- 
ing 0, obtaining j3 = (3.35 ± 0.08) x 10" 3 /r7i 7r+ where 
Re£ + hp -> K+n) = (28.06 ±0.27 ±0.45) x lO^/m^ 
lfl6ll and a (tt+ti tt P ) = - (0.1195 ± 0.0016)x 10 _3 /m w + 



■ 1711 - In HBChPT up to one loop and O (q 4 ), ft is fixed by 
the imaginary part of Eo + — that is parameter- free — providing 
PhbciiPT = 2.71 x 10~ 3 /m 7r + which is far away from the 
unitary value. Because of the lack of unitarity of the S-wave 
amplitude [ 10] it is customary to substitute the S wave provided 
by HBChPT by a unitary prescription Jgl [l^, [111 . However, in 
this Letter instead of substituting the entire S wave for a pre- 
scription we prefer to substitute only the cusp part in Eq + from 
HBChPT by the cusp part of Eq + in equation (|8), keeping the 
smooth part provided by HBChPT. In this way we keep the Eq + 
counter-term and both HBChPT and U-HBChPT approaches 
have the same LECs to fit to the data. 

2.3. Empiricalfit 

The empirical fit is parameterization of the S and P waves 
with a minimal physics input: unitarity in the S wave through 
the (3 parameter and the angular momentum barrier. This is 
accomplished with a Taylor expansion in the pion energy in the 
center of mass u up to first order on the smooth part of Eq + 
and Pi/q adding the cusp part in equation ([H) to the S wave 
and keeping the imaginary part of the P waves equal to zero, in 
summaiyQ 
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where <\ <>, P<°\ p[ 1] , P 2 (0) , P 2 (1) , P 3 (0) , and P3 are 
free parameters that will be fitted to the experimental data. We 
note that this expansion goes to a lesser order in cj than HBChPT 
- i.e. Eq 1 + in equation Q goes to order w 3 - but entails more 
parameters. We note that chiral symmetry is not imposed in this 
approach. 



1 The empirical parameterization in @ 13 expands on the photon energy in 
the laboratory frame E 1 while we prefer to expand in the pion energy in the 
center of mass frame u) in order to have direct comparison to HBChPT. Both 
approaches render equally good description of the observables and provide the 
same multipoles. 
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Figure 1: (Color online.) x 2 /dof energy dependence for the empirical (full 
black squares), HBChPT (full green circles), and U-HBChPT (open blue cir- 
cles) fits. Each point represents a separate fit and the connecting lines are 
drawn to guide the eye. The points are plotted at the central energy of each 
bin, although the calculations take the energy variation inside of each bin into 
account. The value x 2 /dof= 1 is highlighted with a solid line. 



Figure 2: (Color online.) Upper energy (fit) dependence of the LECs. Error 
band and full circles (green): HBChPT fits; Empty circles (blue): U-HBChPT 
fit. The high correlation between LECs a+ and a_ makes their error bars 
larger. Errorbars for the U-HBChPT LECs are not depicted but they are ap- 
proximately of the same size as the HBChPT ones. 



3. Results 

Equipped with the HBChPT, U-HBChPT and empirical ap- 
proaches we perform fits to the experimental data in up 
to different maximum photon energies _E™ ax within the range 
[158.72, 191.94] MeV and compute the x /dot as well as the 
corresponding errorbars of the extracted parameters (see 
The energy bins of the data are approximately 2.4 MeV wide, 
which is taken into account in the fitting and calculations. We 
do not employ the first two energy bins from 0, 146.95 and 
149.35 MeV, because they are less reliable due to systematic 
errors. The amount of data employed in each fit depends on 
up to what energy we are fitting, — i.e. for our lowest-energy 
fit (£™ ax = 158.72 MeV) we employ 100 experimental data 
(80 differential cross sections and 20 photon beam asymme- 
tries) and for our highest-energy fit (_E™ ax = 191.94 MeV) we 
employ 514 experimental data (360 differential cross sections 
and 154 photon beam asymmetries). The highest-energy fit has 
been chosen high enough to obtain a x 2 /dof that ensures that 
the three approaches no-longer hold and the lowest-energy fit to 
ensure a reliable fit with enough experimental data. Systematics 
are not included in the \ 2 and this uncertainty can amount up to 
4% in the differential cross section and 5% in the photon asym- 
metry. The fits are performed employing a genetic algorithm 
whose details can be found in 1 18]. 

3.1. Quality of the fits 

Figure Q] shows the x 2 /dof for every fit performed versus 
the upper energy i?™ ax of the fit as well as the number of data. 
It is shown that up to ~170 MeV all the fits are equally good 
providing very low x 2 /dof. Above 170 MeV the trend is differ- 
ent; while the empirical fit remains with a good and stable x 2 / 
dof, both the HBChPT and the U-HBChPT start rising, a trend 
that shows clearly how the theory fails to reproduce the experi- 
mental data above that energy. Because we obtain very similar 



result for U-HBChPT and HBChPT, lack of unitarity cannot be 
blamed for the disagreement between theory and experiment. 
The HBChPT result contrasts with the empirical fit that up to 
180 MeV provides a good description of the data. Above 185 
MeV the x 2 /dof of the empirical fit starts to rise showing the 
effects of higher orders in the partial waves and the appearance 
a non-negligible contribution from the imaginary part of the 
i waves. 

3.2. LECs as a function of E" ,ax 

An important test of the accuracy of the HBChPT expan- 
sion is the stability of the empirical LECs versus E™ m . The 
empirical fit provides a solid benchmark because the parame- 
ters are the same (within errors) in the whole energy region. 
Figure |2] shows the £"™ ax (fit) dependence of the LECs for both 
the HBChPT (with errors) and U-HBChPT approaches. This 
includes the S-wave LECs a + and a_ in Figures[2(a) and[2(b) 
respectively and P-wave LECs £i, £2, and b p in Figures 0(c), 
|2](d), and|2(e). Errors are larger for the fits with lowest E™ x 
because of the smaller amount of data. The S-wave LECs are 
fairly stable in the whole energy range and both HBChPT and 
U-HBChPT are approximately constant within errors. On the 
contrary, P-wave LECs show a non-stable pattern with a posi- 
tive slope for £1 and b p and a negative slope for £2- The large er- 
ror bars make the extracted LECs compatible up to i?™ ax ~ 175 
MeV except for £ x , whose value for the £™ ax = 170.53 fit is al- 
ready incompatible with the lower energy fit £"™ ax = 161.08, 
confirming that ~170 MeV is the upper energy limit of the 
theory. Besides, approximately at ^170 the U-HBChPT and 
HBChPT P-wave LECs start to be incompatible. The U-HBChPT 
LECs are systematically smaller in absolute value than the ones 
obtained through HBChPT, this is expected because the unitary 
P is larger than Phbckpt giving a larger contribution by the 
Im£'o+ which has to be compensated by the other multipoles. 
The slopes of the P-wave LECs show that higher order, rela- 
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Figure 3: (Color online.) Real part of the S and P waves. HBChPT fit up to 
168.16 MeV: Error band and solid line (green); Empirical fit up to 180.02 MeV: 
Dashed. U-HBChPT lit up to 168.16 MeV: Dash-dotted; The data shown are 
the single energy multipoles extracted from the experimental differential cross 
sections and asymmetries in 01. The gray area above the energy axis repre- 
sents the systematic errors |;4] and the red area at the top of (a) the uncertainty 
associated to our knowledge of the D waves. 



tivistic and A(1232) effects are absorbed into them, calling for 
improvement in the theory. 

We have also looked into the correlations by computing the 
correlation coefficient r(x, y) = a xy / (a x a y ) for each pair of 
parameters and for every HBChPT and U-HBChPT fit. We 
find that the correlation remains more or less stable for each 
pair throughout every fit. The S-wave LECs are highly corre- 
lated r(a + ,a_) « [0.78,0.88], £2 and b p provide r(£ 2 ,& P ) ~ 
[0.45, 0.6], which is not unexpected due to the photon asymme- 
try Ws response structure |£l UH, and the rest are fairly uncor- 
rected lying within the range [—0.2, 0.2]. In the case of the S 
wave the correlation is responsible of the large error bars asso- 
ciated to a+ and et_ and indicates that energy dependence and 
threshold value of Rei?o+ cannot be obtained separately with- 
out further experimental information. Regarding the magnitude 
of the empirical LECs, the empirical values are within a factor 
of two of the values estimated in Refs. iTToL IT ] through reso- 
nance saturation. 
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Figure 4: (Color online.) Differential cross section and photon beam asym- 
metry S. Error band and solid line (green): HBChPT fit up to 168.16 MeV; 
Dashed: empirical fit up to 180.02 MeV. The U-HBChPT calculation is not de- 
picted because it completely overlaps with the HBChPT. For the photon asym- 
metry, (c) and (d), the three approaches are undistinguishable. 



3.3. Comparison with experimental data 

In order to compare with experimental data one has to choose 
a best fit (set of parameters/LECs) for each approach. In our 
case we pick the fit up to 180.02 MeV for the empirical fit (x 2 / 
dof= 1.21) and the fits up to 168.16 MeV for HBChPT ( X 2 / 
dof= 1.25) and U-HBChPT (x 2 /dof= 1.21) approaches. Fig- 
ure|3]shows the single energy multipoles extracted from experi- 
mental data compared to the three approaches. The HBChPT 
fit is shown as an error band. The procedure to obtain the 
single energy multipoles from the data is explained in |4|] and 

The 



the error bars are computed as described in Appendix A 
data below the unitary cusp (146.95 and 149.35 MeV) are not 
reliable enough to accurately extract the single energy multi- 
poles and, therefore, are not shown. Overall, the HBChPT 
and U-HBChPT do a reasonable job describing the multipoles. 
However, the energy dependence of RePi / q deviates somewhat 
from the empirical values. Systematically U-HBChPT P waves 
are smaller in absolute magnitude than those extracted through 
HBChPT. This is a consequence of the different (3 value as ex- 
plained in Section 13.21 The discrepancy at threshold between 
HBChPT and unitary fits for Re£'o+ is also due to the value of 

Figure|4]compares the empirical, HBChPT and U-HBChPT 
approaches to the differential cross section and photon beam 
asymmetry at two different energies, one within the HBChPT 
and U-HBChPT fitting region (£ 7 = 161.08 MeV) and an- 
other outside it {E 1 = 175.22 MeV). Two results are notewor- 
thy. First, the photon asymmetry is well reproduced for both 
energies by all the fits, Figures |4](c) and|4](d); if we compare 
with other energies -higher, lower and intermediate- we find 
the same level of agreement between theory and data, obtain- 
ing that all the approaches are of the same quality and provide 
a good description of the photon asymmetry in the whole en- 
ergy range considered in this Letter. Second, the HBChPT and 
U-HBChPT approaches underestimate the cross section for en- 
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Figure 5: (Color online.) Components of the Wt and Ws responses for the 
HBChPT fit up to 168. 16 MeV (green error band and solid line) and empirical 
fit up to 180.02 MeV (dashed). All the errors are computed at a x 2 + 1 level. 
If the error bars are increased to 2<r the Si extraction is compatible with zero 
in the whole energy range. 



ergies above the fitting limit (i?™ x = 168.16 MeV), as can be 
seen in Figure|4](b). This situation is clearer if we look into the 
component To of the differential cross section response Wt in 
Figure[5](a) which above ^170 MeV is largely underestimated 
by the HBChPT approach. The significant components of the 
Wt response To, T\, and T2 are obtained fitting the differential 
cross sections for each energy bin to equation ([1}. The same 
operation is done with the photon asymmetry, fitting the data 
to equation (O, extracting So and Si. The other two signifi- 
cant components of the Wt response T\ and T2 are fairly well 
described up to ~ 175 MeV. In the case of Ws, So is well deter- 
mined and HBChPT describes it fairly well in the whole energy 
range. 

3.4. Probing D waves 

The Si component of Ws in Figure[5](e) is due to the inter- 
ference among P and D waves JHEa]. The empirical values of 
Si are consistent with the Born terms contribution of D waves 
but unfortunately the experiment is not accurate enough in the 
photon asymmetry to provide a quantitative measurement. The 
small non-zero effect between 175 and 185 MeV disappears if 
errors are computed at a 2a level. Hence, with the current ex- 
perimental information only four quantities can be accurately 
obtained from each energy bin To, T\, and T2 from the differ- 
ential cross section and So from the photon asymmetry, that 
allows to extract the four multipoles in Figure [3] If we intend 
to obtain information on the rest of the multipoles we need ei- 
ther more accuracy in the data to pin down Si (to obtain infor- 
mation on D waves), T3 (PxD interference) or T4 (D waves), 
or to measure other observables like the target asymmetry (to 
obtain lmEo + ) Il9l 2011. the E asymmetry ||9|] (D waves) or 
the F asymmetry I® [l9i |20|l (ReE 0+ and D waves). With our 
current knowledge of the P waves, the accurate extraction of 
Imi?o+ and the (3 parameter from the target asymmetry is feasi- 
ble. This observable has been measured at Mainz together with 



the F asymmetry and data analysis is currently in the progress 
Ml 

Returning to D waves, they have been incorporated in our 
analysis as the Born terms and the Si component of Ws is con- 
sistent with this approach. Up to order 0(q 4 ) in HBChPT this 
is the only contribution together with an E2- counter-term i2lll 
which provides an additional LEC. However we have neglected 
it in our calculation because it has no impact in the x 2 , and, 
therefore, it cannot be determined. Future more accurate data 
will provide a measure of the D-wave effects and allow to pin 
down the E2- counter-term if a deviation from Born terms is 
found. 

4. Conclusions 

Because of the high-quality experimental data gathered by 
the A2 and CB-TAPS collaborations at Mainz we can asses the 
electromagnetic multipoles and their energy dependence to the 
best precision ever and we can accurately test chiral symmetry 
and the range of application of approaches like HBChPT. Based 
on the accumulated evidence -LECs stability, x 2 /dof and the 
empirical fit which works up to ~180 MeV- we can establish 
a clear upper-energy limit of validity for HBChPT. We have 
found this limit to be ~ 170 MeV of photon energy in the labo- 
ratory frame. The slopes of the P-wave LECs in Figure |2]show 
that higher order, relativistic and A(1232) effects are absorbed 
into them, calling for improvement in the theory. The lack of 
unitarity in the S wave is not responsible of the disagreement 
between HBChPT and the experimental data as we have proved 
through the U-HBChPT approach. Besides, we do not think 
that calculating higher orders in HBChPT will help extending 
the energy range of application of the theory because data are 
fairly well reproduced by the empirical fit up to 180 MeV which 
expands up to a lower order in pion energy. Further improve- 
ment in the theory is necessary if chiral perturbation theory is to 
be applied above 170 MeV. However, the relativistic calculation 
does not provide better agreement with data than the HBChPT 
approach fl l^ll]. T his suggests the inclusion of the A(1232) in 
the calculation H22I1 . 

The fact that we have achieved accurate knowledge on the 
P-wave multipoles allows us to undertake the study of other 
multipoles (D waves and Im£'o+) and observables (target and F 
asymmetries) through further experimental research that is un- 
derway ll20ll . Current experimental information does not allow 
to test our knowledge on D waves but we are hopeful about 
forthcoming experiments. 

Acknowledgements 

We thank the A2 and CB-TAPS collaborations for mak- 
ing available the experimental data prior to publication. C.F- 
R. is supported by "Juan de la Cierva" programme of Spanish 
Ministry of Economy and Competitiveness and his research has 
been conducted with support by Spanish Ministry of Economy 
and Competitiveness grant FIS2009-1 1621-C02-01, the Mon- 
cloa Campus of International Excellence (CEI Moncloa), and 



5 



by CPAN, CSPD-2007-00042Ingenio2010. A.M.B. research is 
supported in part by the US Department of Energy under con- 
tract No. DE-FC02-94ER408 18. 

Appendix A. Error bar calculation 

Error bars have been computed through a Monte Carlo (MC) 
simulation. Once the minimum Xmin nas been assessed the x\ 
defined as xi — Xmin + 1- Once we have the x\ we run a 
MC varying the values of the parameters, we compute the cor- 
responding x 2 f° r eacn set °f parameter values, and we accept 
those sets which provide x 2 < Xi- If enough statistics are col- 
lected, the boundary of the simulation defines the confidence el- 
lipse and the error bars for each parameter lEil . We also obtain 
correlation plots between parameters as well as the correlation 
matrix. Once the MC has been run we have a file with thou- 
sands of combinations of the parameters which are within the 
Xi level. We use those sets to compute the bands in the partial 
waves and the observables which are shown in the figures. In 
this way the error bands in the partial waves and the observables 
take properly into account the correlations among parameters. 



[22] T.R. Hemmert, B.R. Holstein, J. Kambor, J. Phys. G 24 (1998) 1831; 
Phys. Lett. B 395 (1997) 89; V. Lensky, V. Pascalutsa, Eur. Phys. J. C 65 
(2010) 195; J. M. Alarcon, J. Martin Camalich, J.A. Oiler, Phys. Rev. D 
85 (2012) 051503; J. M. Alarcon, J. Martin Camalich, J.A. Oiler, submit- 
ted for publication (2012), arXiv: 1210.4450 [hep-ph]. 

[23] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86 (2012) 010001; 
G. Cowan, Statistical Data Analysis, Oxford University Press, Oxford 
2002. 



References 

[1] J. F. Donoghue, E. Golowich, B. R. Holstein, Cambridge Monographs in 
Particle Physics, Nuclear Physics and Cosmology Vol. 2: Dynamics of 
the Standard Model, Cambridge University Press, Cambridge 1992. 

[2] J. Gasser, M. E. Sainio, A. Svarc, Nucl. Phys. B 307 (1988) 779. 

[3] V. Bernard, N. Kaiser, U.-G. MeiBner, Int. J. Mod. Phys. E 4 (1995) 193. 

[4] D. Hornidge et al, submitted for publication (2012), arXiv : 121 1.54951 
[nucl-ex]. 

[5] V. Bernard, N. Kaiser, J. Gasser, U.-G. MeiBner, Phys. Lett. B 268 (1991) 
291; V. Bernard, N. Kaiser, U.-G. MeiBner, Nucl. Phys. B 383 (1992) 
442; V. Bernard, U.-G. MeiBner, Annu. Rev. Nucl. Part. Sci., 57 (2007) 
33. 

[6] A. M. Bernstein, S. Stave, Few Body Syst. 41 (2007) 83. 

[7] A. M. Bernstein, E. Shuster, R. Beck, M. Fuchs, B. Krusche, H. Merkel, 

H. Stroher, Phys. Rev. C 55 (1997) 1509. 
[8] C. Fernandez-Ramirez, A. M. Bernstein, T. W. Donnelly, Phys. Lett. B 

679 (2009) 41; C. Fernandez-Ramirez, PoS(CD09) (2009) 055. 
[9] C. Fernandez-Ramirez, A. M. Bernstein, T. W. Donnelly, Phys. Rev. C 80 

(2009)065201. 

[10] V. Bernard, N. Kaiser, U.-G. MeiBner, Z. Phys. C 70 (1996) 483. 

[11] V. Bernard, N. Kaiser, U.-G. MeiBner, Eur. Phys. J. A 11 (2001) 209. 

[12] P. Buttiker, U.-G. MeiBner, Nucl. Phys. A 668 (2000) 97. 

[13] A. M. Bernstein, Phys. Lett. B 442 (1998) 20. 

[14] B. Ananthanarayan, Phys. Lett. B 634 (2006) 391. 

[15] M. Hoferichter, B. Kubis, U-G. MeiBner, Phys. Lett. B 678 (2009) 65. 

[16] E. Korkmaz et al, Phys. Rev. Lett. 83 (1999) 3609. 

[17] V. Baru, C. Hanhart, M. Hoferichter, B. Kubis, A. Nogga, D.R. Phillips, 

Nucl. Phys. A 872 (201 1) 69. 
[18] C. Fernandez-Ramirez, E. Moya de Guerra, A. Udias, J. M. Udias, Phys. 

Rev. C 77 (2008) 065212. 
[19] A. M. Bernstein, M. W. Ahmed, S. Stave, Y. K. Wu, H. R. Weller, Annu. 

Rev. Nucl. Part. Sci. 59 (2009) 115. 
[20] D. Hornidge, A.M. Bernsein, Eur. Phys. J.: Special Topics 198 (2011) 

133; S. Schumann, AIP Conf. Proc. 1441 (2012) 287; A. M. Bernstein et 

al, Mainz Exp. A2/10-2009, Measurement of Polarized Target and Beam 

Asymmetries in Pion-Prodution on the Proton: Test of Chiral Dynamics 

(2009). 

[21] M. Hilt, Photo- and Electro-Pion Production in Chiral Effective Field 
Theory, PhD dissertation, University of Mainz (201 1); M. Hilt, S Scherer, 
L. Tiator, manuscript in preparation (2012). 



6 



